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ABSTRACT 

Wc present an analysis of the Hubble diagram for 12 normal type la supernovae (SNe 
la) observed in the near-infrared J and H bands. We select SNe exclusively from the 
redshift range 0.03 < z < 0.09 to reduce uncertainties coming from peculiar velocities 
while remaining in a cosmologically well-understood region. All of the SNe in our sam- 
ple exhibit no spectral or S-band light-curve peculiarities and lie in the S-band stretch 
range of 0.8 - 1.15. Our results suggest that SNe la observed in the near-infrared (NIR) 
arc the best known standard candles. We fit previously determined NIR light-curve 
templates to new high-precision data to derive peak magnitudes and to determine the 
scatter about the Hubble line. Photometry of the 12 SNe is presented in the natural 
system. Using a standard cosmology of (iJo, f^m, ^^a) ~ (70,0.27,0.73) we find a me- 
dian J-band absolute magnitude of Mj — —18.39 with a scatter of crj = 0.116 and a 
median if-band absolute magnitude of Mh = —18.36 with a scatter of cth — 0.085. 
The scatter in the H band is the smallest yet measured. We search for correlations 
between residuals in the J- and iJ-band Hubble diagrams and SN properties, such 
as SN colour, B-band stretch and the projected distance from host-galaxy centre. 
The only significant correlation is between the J-band Hiibblc residual and the J-H 
pseudo-colour. We also examine how the scatter changes when fewer points in the 
near-infrared are used to constrain the light curve. With a single point in the H band 
taken anywhere from 10 days before to 15 days after B-band maximiim light and a 
prior on the date of iJ-band maximum set from the date of B-band maxiirnim, we 
find that we can measure distances to an accuracy of 6%. The precision of SNe la in 
the NIR provides new opportunities for precision measurements of both the expansion 
history of the universe and peculiar velocities of nearby galaxies. 

Key words: cosmology: observations, (cosmology:) distance scale. 
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1 INTRODUCTION 

SNe la are widely accepted to be excellent standardizable 
candles at optical wavelengths. The development of various 
empirical corrections, which can relate peak luminosity and 



sensitivity variations. We used the XDIMSUM package in 
IRAF to remove the sky and a modified version of ISIS2 



light-curve shape ( 


Phillips 1993 Hamuy et al. 19961, SN 


colour (Tripp 1998 


Jha et al. 2007 


1, spectral information 


(Bronder et al. 2007 Walker et al. 


2011 Silverman et al. 


20121 and host-galaxy mass (Sullivan et al. 20101, allow 



SNe la peak magnitudes to serve as distance indicators 
with a corrected scat teiP] of as low as 0.13 magnitudes 
( Silverman et al.||2012[ |Conley et al.|2011 |. 

There has been growing evidence, though, that SNe 
la may be more accurate in the near-infrared (NIR) (for 
a recent review, see Phillips 20111. Near-infrared light 



experiences less attenuation from dust than light at optical 
wavelengths, and theoretical models predict a smaller 
intrinsic dispersion in the near-infrared peak magniude 
( Kasen |2006j . Furthermore, there appears to be little or no 



relationship between light-curve shape and peak luminosity 
( Wood-Vasey et al.|2008[|Folatelh et al.|2010|pandel et al 



2011 1, meaning no empirical corrections need be applied. 



Recent studies (Meikle 


2000 


Wood-Vasey et al. 


2008 


FolateUi et al.||2010| Krisciunas et al.||2004| |Kattner et al.| 


2012 Mandel et al. 


2011 


have shown that SNe la in 



the NIR can be as reliable as corrected SNe la observed 
at optical wavelengths, but these studies have had some 
limitations. Foremost among these limitations is that the 
SNe were observed at distances that are not sufficiently 
large to place them in the Hubble flow (defined here as 
z > 0.03) and hence are affected by peculiar velocities. Here 
we present a study of a sample of Type la SNe selected 
to lie exclusively in the redshift range 0.03 < z < 0.09 
in order to minimise redshift uncertainty due to peculiar 
velocity. We use NIR images on 8m class telescopes to 
observe supernovae in our sample, resulting in photometry 
that is more precise than that achieved for the handful of 
supernovae that have been observed beyond z = 0.03. 

2 SNE lA SAMPLE 

The SNe la used in our sample are summarised in Table 
[2] These SNe la were discovered by the Palomar Transient 
Factory (PTF; |Rau et al.||2009[ [Law et al.||2009[), spectro- 



scopically confirmed, observed in the p, r, i bands (Maguire 
et al. in prep.) at the Liverpool Telescope and followed up 
in both the NIR J and H bands, each with four epochs 
of observation. PTF09dlc was observed using HAWK-I 
( |CasaU et al.|2006| l on ESO's 8.1m VLT, and all others were 
observed using the NIR Imager and Spectrometer 

( Hodapp et al.|2000[ ) on Gemini Observatory's 8.2m Gemini 
North. Reference images were taken ~ 1 year after initial 
observations so that host-galaxy light could be removed. 

The data were processed in a standard manner using 
IRAF0 and our own scripts. Darks were used to remove 
the pedestal, and flats were used to remove pixel-to-pixel 



^ We use the rms about the mean as a measure of the scatter. 
^ http:/ /www. eso.org/sci/facilities/paranal/instruments/hawki/ 
and http://www.gemini.edu/sciops/instrumcnts/niri/ 
^ IRAl"' is distributed by the National Optical Astronomy Ob- 
servatories which are operated by the Association of Universities 



( Alard 2000 1 to subtract the reference image of each SN 
from the images with SN light, leaving just light of the 
SN. Zero-points were derived from standards found in the 



Persson et al. ( 1998 1 catalogue. All magnitudes are reported 
in the natural system of each instrument. K-corrections 
were calculated using the revised spectral template of |Hsiao| 
et al. (20071 and corrections for galactic dust extinction 



were applied. We did not warp the spectral template to 
match the observed colour. 

For the SNe PTF09dlc, PTFlOhmv, PTFlOmwb, 
PTFlOnlg, PTFlOtce, PTFlOufj, PTFlOwnm and 
PTFlOxyt, the first observing epoch was before NIR 
maximum. For the remaining SNe, the first observing epoch 
was near NIR maximum. When fitting light-curves for the 
SNe, we included the date of B-band maximum as an extra 
prior. The i?-band maximum was determined using SiFTO 



(Conley et al. 20081, as was the B-band stretch (Maguire 
et al. in prep.). 

The _B-band stretches range from 0.8 to 1.15, so there 
are no fast decliners in our sample. Nor are there any 
heavily reddened supernovae. Furthermore, all of the SNe 
appear to be spectroscopically normal. The absence of fast 
decliners and heavily reddened SNe is a potential limitation 
of the current sample. 

The photometry for all SNe is presented in Table [l] 

3 TEMPLATE FITTING 

To find the light-curve maxima, we flt the Flexible Light- 



curve InfraR ed Template (FLIRT) presented in Mandel 
et al. ( 2009 1 to our K-corrected photometrjj^ As we did 
not observe all SNe la before maximum light, we fit the 
template to the data with an additional constraint. Concur- 
rent observations of SNe la at NIR and optical wavelengths 
provide evidence that NIR maxima occur ~ 5 days before 
_B-band maximum ( |Meikle|2000l [Krisciunas et al.|2004| ). We 
use this prior assumption as a constraint on the date of NIR 
maximum, tp, with a 1—a error of cTp = 1 day, imposed on 
the likelihood function. By fitting the light-curve template 
to the SNe la observed before maximum without a prior, 
we determine that the J-band maximum occurs 5.36 days 
before the _B-band maximum with a standard deviation 
of 0.74 days and the H-hand maximum occurs 4.28 days 
before with a standard deviation of 0.70 days. Figure [2] 
displays the template-fitted light-curve for PTFlOtce. 

The errors in the peak apparent magnitudes in both 
the J and H band were calculated via 100 Monte Carlo 
simulations of template fitting to the data, including 
repeated Gaussian-distributed random samplings within 
the errors of each point. 

To convert apparent magnitudes to absolute mag- 
nitudes, we determine the distance modulus for each 
SN, calculated using the standard fiat cosmology of 
(Ho,Slm,r2A) = (70,0.27,0.73). The redshift measurement 
for PTFlOufj is uncertain since the redshift is measured 
from the supernova and not the host. In this case, we have 



for Research in Astronomy, Inc. under the cooperative agreement 
with the National Science Foundation 

* Photometry to be published in Barone-Nugent et al. (in prep.). 
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NIR type la SNe 3 



SN name 


MJD 


mj 










PTFOQdlc 


55068 


14 


19.803 ± 0.023 


18 


813 







021 




55073 


11 


19.906 ± 0.040 


18 


933 


± 





044 




55085 


16 


20.997 ± 0.059 


19 


646 


± 





050 




55093 


18 


21.399 ± 0.067 


19 


635 


± 





028 


PTFlOhdv 


55340 


25 


19.084 ± 0.044 


18 


516 


± 





052 




55346 


25 


19.407 ± 0.057 


18 


766 







052 




55350 


25 


19.787 ± 0.085 


18 


920 







065 




55354 


41 


20.266 ± 0.139 


19 


048 


± 





085 


PTFlOhmv 


55341 


26 


- 


17 


887 


± 





019 




55346 


28 


- 


17 


501 


± 





026 




55351 


30 


- 


17 


595 


± 





017 




55357 


25 


- 


17 


823 


± 





054 




55358 


25 


- 


17 


825 


± 





019 




55359 


26 


- 


17 


826 







018 


PTFlOmwb 


55378 


41 


19.066 ± 0.031 


18 


429 







030 




55381 


43 


18.677 ± 0.027 


17 


765 







014 




55386 


41 


17.683 ± 0.026 


17 


484 


± 





012 




55391 


43 


18.202 ± 0.022 


17 


569 







013 




55396 


39 


18.584 ± 0.023 


17 


760 







014 


PTFlOndc 


55386 


43 


19.898 ± 0.062 


19 


249 







079 




55389 


38 


20.110 ± 0.088 


19 


468 


± 





099 




55393 


39 


20.333 ± 0.092 


19 


488 


± 





107 




55395 


35 


20.496 ± 0.101 


19 


519 


± 





118 


PTFlOnlg 


55386 


39 


19.352 ± 0.044 


18 


694 







039 




55387 


33 


- 


18 


500 


± 





041 




55389 


35 


- 


18 


571 


± 





035 




55393 


36 


19.679 ± 0.073 


18 


770 


± 





044 




55395 


37 


19.816 ± 0.059 


18 


791 


± 





057 


PTFlOqyx 


55422 


63 


19.652 ± 0.069 


19 


005 


± 





046 




55424 


56 


19.701 ± 0.045 


19 


121 


± 





049 




55428 


64 


20.156 ± 0.076 


19 


487 







102 




55433 


57 


20.654 ± 0.104 


19 


736 







088 


PTFlOtce 


55435 


41 


18.634 ± 0.040 


18 


047 


± 





032 




55438 


34 


18.592 ± 0.036 


17 


943 


± 





026 




55443 


50 


18.778 ± 0.041 


18 


164 







028 




55448 


53 


19.175 ± 0.059 


18 


352 







035 


PTFlOufj 


55447 


55 


20.496 ± 0.112 


19 


444 


± 





079 




55450 


64 


20.387 ± 0.093 


19 


241 


± 





068 




55455 


50 


20.312 ± 0.100 


19 


202 


± 





061 




55460 


58 


20.996 ± 0.238 


19 


544 


± 





097 


PTFlOwnm 


55469 


38 


19.513 ± 0.055 


18 


952 


± 





044 




55471 


48 


19.443 ± 0.052 


18 


822 


± 





040 




55476 


27 


19.482 ± 0.053 


18 


984 


± 





069 




55481 


28 


19.883 ± 0.072 


19 


262 


± 





126 


PTFlOwof 


55470 


33 


19.086 ± 0.030 


18 


488 


± 





033 




55473 


29 


19.201 ± 0.039 


18 


592 


± 





042 




55478 


24 


19.491 ± 0.063 


18 


855 


± 





051 




55483 


29 


20.071 ± 0.071 


19 


173 


± 





061 


PTFlOxyt 


55486 


26 


19.116 ± 0.051 


19 


071 


± 





118 




55492 


24 


19.277 ± 0.047 


19 


438 


± 





179 




55496 


21 




19 


735 







238 




55498 


20 




19 


596 


± 





227 



Table 1. The J- and //-band photometry for each epoch of each SNe in our sample. Apparent magnitudes have been corrected for 
Galactic dust extinction but not K-corrected. 



assumed an uncertainty of = 0.005 (Maguire et al. in 
prep.). For all other redshifts, the errors make negligible 
difference to the distance modulus. 

4 RESULTS 

The peak apparent magnitudes for each SN la in the J and 
H band are quoted in Table [2] In the H band, the me- 
dian absolute magnitude of the supernovae in our sample is 



Mh ~ — 18.36±0.04. This compares reasonably well with the 



median magnitudes of supernovae in Folatelli et al. (20101, 



Wood-Vasey et al. ( 2008 1 and Krisciunas et al. ( 2004 1 for 



which we derive Mh = -18.37 ± 0.07, Mh = -18.29 ± 0.04, 
and Mh = —18.43 ± 0.06, respectively. The errors on the 
median are computed by jackknife resampling. In part, the 
difference in magnitudes comes from the limited number of 
supernovae in each of the samples. The differences may also 



© 0000 RAS, MNRAS 000, 000-000 



4 Barone- Nugent et al. 



Supernova 


RA(J2000) 


DEC (J2000) 


^helio 


M B 


max 




PTF09dlc 


21:46:35.5 


+06:23:23.5 





0678 


0.0666 


i9.oio±i;-»i« 


lO. iOZ Q Q2Q 


PTFlOhdv 


12:07:43.4 


+41:29:27.9 





0534 


0.0542 


1« K91 +0.031 
03T 


18.453+°-°26 

U.U.iD 


PTFlOhmv 


12:11:33.0 


+47:16:29.8 





0320 


0.0327 


17.431+^0^ 




PTFlOmwb 


17:17:50.0 


+40:52:52.1 





0313 


0.0312 


„ +0.007 
± 1 .ot J^_0.007 


17 on+0 010 
± t .^i-0_Q 009 


PTFlOndc 


17:19:50.2 


+28:31:57.5 





0818 


0.0817 


19.330+0511 


1Q 007+0.051 
I3.00/_o 


PTFlOnlg 


16:50:34.5 


+60:16:35.0 





0560 


0.0559 


18.634lH3,« 


18.720in29 


PTFlOqyx 


02:27:12.1 


-04:31:04.8 





0660 


0.0652 


19.122lH^^ 


19.016irof3 


PTFlOtce 


23:19:11.0 


+09:11:54.2 





0410 


0.0398 


17.988lHi? 


17.868iro^? 


PTFlOufj 


02:25:39.1 


+24:45:53.2 





077 ± 0.005 


0.0762 ± 0.005 


19 O7fi+0 036 


19.306in^? 


PTFlOwnm 


0:13:47.26 


+27:02:26.0 





0656 


0.0645 


18.94llHf7 


18.764irog 


PTFlOwof 


23:32:41.8 


+15:21:31.7 





0526 


0.0514 




18.458iroS 


PTFlOxyt 


23:19:02.4 


+13:47:26.8 





0496 


0.0484 


18.379tHf3 


i8.45oiror7 



Table 2. The sample of SNe la showing their redshifts and J- and J?-band peak apparent magnitudes. PTFlOndc, PTFlOtce, PTFlOufj, 
PTFlOwof, PTFlOxyt were located by the citizen science Galaxy Zoo: Supernovae project ( [Smith et al.|20TT| |. 




Figure 1. The upper panels shows the Hubble diagrams of our sample of 12 SNe (green circles) in the H band (left) and the J band 
(right), including the |Wood-Vasey et al.| l |2008[ l, [Folatelli et al.| | [2010[ l and |Krisciunas et al.| | [2"004[ | samples (blue crosses) for comparison. 
The red line represents the apparent magnitude that would be observed assuming a constant absolute magnitude of —18.36 (H band) 
and —18.39 (J band). The lower panels show the Hubble residual, i.e. the deviation from the red line. The shaded red region, z < 0.03, 
is the region excluded in our sample due to the associated peculiar velocity errors. The solid black lines represent the change in the 
distance modulus due to a peculiar velocity of ±300 kms~^ 



be due to the difTerences in the photometric systems that 
each author adopts. These differences are generally smaller 
than 0.04 mag. 

In the J band, the variations between samples are larger 
and more difficult to understand. We find a median magni- 
tude of Mj — —18.39 ± 0.06, compared to median mag- 
nitudes of Mj = -18.48 ± 0.05, Mj = -18.2 9 ± 0.10, 
and Mj = —18.73 ± 0.05 for the supern ovae in |Folatelli 



(20041. The errors were calculated as for the H band. The 



et al. (20101, Wood-Vasey et al. (20081 and Krisciunas et al 



largest source of uncertainty is due to sample size. Difference 
due to photometric systems are generally smaller. In order 
to constrain the dark energy equation of state beyond limits 
currently derived using supernovae observed in the optical, 
the mean magnitude needs to be determined to a level of 
accuracy that is better than 0.01 mag. 

The J- and iif-band Hubble diagrams are shown in 
Figure [T] Our results are compared with the surveys of 
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Figure 2. The FLIRT template fitted to the light-curves for PTFlOtce in the H band (left) and the Jband (right). 



Wood-Vasey et al. (20081 and Folatelli et al. (20101. Our 


X 








0"H,tot 


observed scatter about the mean absolute magnitude of 


1 


0.002 


0.146 


0.007 


0.116 


o-H = 0.085 ± 0.016 mag and aj = 0.116 ± 0.027 mag rep- 


2 


0.012 


0.126 


0.002 


0.096 


resent a considerable improvement on previous samples, in- 


3 


0.008 


0.120 


0.002 


0.089 


cluding Wood- Vasey et al., (j2008|) who observed cth = 0.28 


4 


0.000 


0.116 


0.000 


0.085 



mag and aj = 0.41 mag, [Krisciunas et al.| ( |20( 

served cth = 0.17 mag and aj = 0.13 mag, [Folatelli et aTl 



(20101 who observed cth = 0.33 mag and ctj = 0.37 mag 
and [Mandel et aT] ( |2009[ ) who observed cth = 0.11 mag and 
Ctj = 0.17 mag. The values for ctj and cth presented here are 
the scatter in the magnitude at J-max and iJ-max, respec- 
tively. 

Photometric errors make little contribution to the scat- 
ter we measure. Peculiar velocities only make a significant 
contribution if they are large (i.e. ~ 300kms~^). The intrin- 
sic scatters we measure, assuming peculiar velocities of 300 
kms~^ are CTint = 0.101 in the J band and CTint = 0.065 in the 
H band. For peculiar velocities of 150 kms~^ we obtain in- 
trinsic scatters of CTint = 0.109 in the J band and CTint ~ 0.077 
in the H band. Part of the scatter we observe may also be 
due to spectral variability aflecting the K-corrections. 

We note that the light-curves of some of the SN in our 
sample are not adequately described by the light-curve tem- 
plate that we have used. This is because there is a diversity 
in light-curve shapes, which may lead to opportunities to 
reduce the scatter about the Hubble line further. We inves- 
tigated this by recomputing the scatter after removing SNe 
with poor light-curve fits. If we eliminate light-curves with a 
X'^ /Degrees of Freedom greater than 5 (a chance probability 
of Q = 0.002), we reject PTFOQdlc and PTFlOmwb in the 
J band and PF09dlc, PTFlOhmv and PTFlOmwb in the H 
band. The resulting values of scatter in the J and H bands 
are 0.105 and 0.097 respectively. The scatter in the H band 
increases marginally but is noticeably lower for the J band. 
PTF09dlc, PTFlOhmv and PTFlOmwb exhibit no spectral 
or B-band light-curve peculiarities. However, applying such 
cuts for the purpose of reducing the scatter needs to be done 
with care since they may introduce biases that may affect 
the accuracy with which cosmological parameters can be 
measured. 

The relationship between the peak magnitude in both 
filters and pseudo-colour (mj^niax-TiH.max) was also consid- 
ered. In the optical, there is a strong correlation between 
colour and luminosity (Tripp 1998 Jha et al. 20071. How- 
ever, we find no obvious correlation between the ff-band 



Table 3. The mean peak magnitude and dispersion in the J and 
H bands when using light-curves with x points. 

residual and Jmax-Hmax pseudo-colour in our sample. On 
the other hand, we do find a correlation between the J-band 
residual and pseudo-colour. From the least-squares residu- 
als in the {/-direction, we find a relation for the line of best 
fit, where the errors on the coefficients are 1 — ct dispersions 
calculated via 1000 Monte Carlo simulations by randomly 
varying the points within their error bars: 



J-band residual = (0.91 ± 0.2) x (i\/j,niax - Mn.max) 

-I- (0.05 ±0.02) 



(1) 



with a scatter about the line of best fit of ct = 0.071. 

We also considered the relationship between the peak 
magnitude in both filters with the angular distance between 
the SN and the centre of the host galaxy and _B-band stretch. 
However, we find no correlations between any of these quan- 
tities. It should be noted that the B-band stretch range cov- 
ered by this sample is limited. 



4.1 Template fitting with a subset of points 

The previous discussion refers to SNe la light-curves con- 
sisting of 4 data points, which we used to fit the light-curve 
templates. In this section we discuss the optimal number of 
epochs that should be observed in the NIR given finite tele- 
scope time. 

We have already calculated the peak magnitude, Mj^h, 
for each SNe la using all data available. We now calculate 
the peak magnitudes using a subset of the data points of size 
X, denoted by M^-pts, and find the scatter of these magni- 
tudes around Mj^h. The total scatter in the sample using x 
points is denoted by ctj, h, tot. The results for these quantities 
based on our SNe la sample are presented in Table |3] 

We find in all cases that there is approximately one- 
hundredth of a magnitude difference in the mean peak mag- 
nitude or less between different numbers of epochs used. 
However there is additional scatter around the mean which 
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X 


Sample size 




SEMh 




SEMj 


1 


12 


0.116 


0.033 


0.146 


0.042 


2 


6 


0.096 


0.039 


0.126 


0.051 


3 


4 


0.089 


0.045 


0.120 


0.060 


4 


3 


0.085 


0.049 


0.116 


0.067 



Table 4. The standard error of the mean for an example of 12 
observing nights distributed among 3, 4, 6 and 12 SNe. 

decreases as more points are added. This illustrates the ex- 
change between amount of telescope time used and disper- 
sion measured. Table [3] suggests that the increased accuracy 
from 3 points to 4 is ~ 0.01 mag or less in both the J and H 
bands. Furthermore, using just one point with a known date 
of maximum can be used while only sacrificing ~ 0.04 mag 
of accuracy (i.e. l-a scatter) in both the J and H bands. 
Importantly, using just one point in the J band gives a dis- 
persion that is similar to the best standardized SNe la ob- 
served in the optical, and using just one point in the H band 
can further improve the accuracy of SNe la observed in the 
optical. We also find that, within the range of epochs cov- 
ered by our data (approximately 10 days before to 15 days 
after _B-band maximum), the amount of additional scatter 
introduced due to using only one point has no dependence 
on when this point is taken with respect to B-band max- 
imum. Importantly, single observations after B-band peak 
are equally as effective as points prior to the peak. 

We may now compare the standard error of the mean 
(SEM) for each subset of points. The SEM for a sample of 
n points is given by. 



SEM = 




where a is the standard deviation of the sample. As an illus- 
trative example, we compare observing 3 SNe over 4 epochs, 
observing 4 SNe over 3 epochs, observing 6 SNe over 2 
epochs and observing 12 SNe for a single epoch. In all cases, 
we assume that we have well-sampled optical light-curves 
for each SN la. Each of these would require equal allocation 
of telescope time. The relative SEM for each method are 
summarised in Table |4] We find that observing four times 
the number of SNe with only a single epoch delivers a lower 
SEM than observing fewer SNe over more than one epoch. 
The larger sample size more than compensates for the higher 
uncertainty in the peak magnitude. 

5 CONCLUSION 

We have presented NIR light curves of 12 SNe la that are 
in the Hubble fiow. We find that the intrinsic scatter in 
peak luminosities of type la supernovae in the NIR J and 
H bands are smaller than previously thought. This is the 
first sample to display an //-band rms scatter as small 
as (Jh = 0.085 ± 0.016 (with a median peak magnitude 
of Mh = —18.36). Our observed J-band rms scatter of 
crj = 0.116 ± 0.027 (with a median peak magnitude of 
Mj = —18.39) is smaller than reported elsewhere. These 
results provide distance errors of ~ 4% using H band 
SNe, making them the most precise standard candles for 
cosmology. 

We have also shown that if concurrent optical obser- 
vations are made, we may use a predicted date of NIR 
maximum as a constraint when fitting the light-curve. With 



this constraint, we may use as few as one NIR observation 
within ~ 5 — 10 days of NIR max per supernova while 
still achieving scatters of cth = 0.116 and oj = 0.146. 
As surveys improve over the coming years and more 
SNe la are discovered, single-night near-infrared observa- 
tions undertaken concurrently with optical observations will 
be the most efficient and accurate way to construct samples. 
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